Cross-sections of spallation residues produced in 1-A GeV 208 Pb on proton reactions 
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Spallation residues produced in 1 GeV per nucleon 208 Pb on proton reactions have been studied using the 
FRagment Separator facility at GSI. Isotopic production cross-sections of elements from 6iPm to g2Pb 
have been measured down to 0.1 mb with a high precision. The recoil kinetic energies of the produced 
fragments were also determined. The obtained cross-sections agree with most of the few existing gamma- 
spectroscopy data. The data are compared with different intranuclear-cascade and evaporation-fission 
models. Drastic deviations were found for a standard code used in technical applications. 

PACS numbers: 25.40.Sc, 24.10.-i, 25.70.Mn, 29.25.Dz 



Spallation reactions have recently captured an increas- 
ing interest due to their technical applications as intense 
neutron sources for accelerator-driven sub-critical reac- 
tors jij or spallation neutron sources g| . The design of an 
accelerator-driven system (ADS) requires precise knowl- 
edge of nuclide production cross-sections in order to be 
able to predict the amount of radioactive isotopes pro- 
duced inside the spallation target. Indeed, short-lived 
isotopes may be responsible for maintenance problems 
and long-lived ones will increase the long term radiotox- 
icity of the system. Recoil kinetic energies of the frag- 
ments are important for studies of radiation damages in 
the structure materials or in the case of a solid target. 
Data concerning lead are particularly important since in 
most of the ADS concepts actually discussed, lead or 
lead-bismuth alloy is considered as the preferred mate- 
rial of the spallation-target. 

The present experiment, using inverse kinematics is 
able to supply the identification of all the isotopes pro- 
duced in spallation reactions and information on their 
recoil velocity. Moreover, the data represent a crucial 
benchmark for the existing spallation models used in the 
ADS technology. The precision of these models to es- 
timate residue production cross-sections is still far from 
the performance required for technical applications, as 
it was shown in Ref. ||. This can be mostly ascribed to 
the lack of complete distributions of all produced isotopes 
to constrain the models. The available data were gener- 
ally obtained by chemistry or gamma spectroscopy 
which give access mostly to cumulative yields produced 
after long chains of decaying isotopes. 

In this Letter we report on complete isotopical pro- 
duction cross-sections for heavy fragments produced in 
spallation of 208 Pb on proton at 1-A GeV, down to 0.1 
mb with a high precision. The kinematic properties of 



the residues are also studied. The cross-sections of lighter 
isotopes produced by fission will be presented in a forth- 
coming publication. 

The experimental method and the analysis procedure 
have been developed and applied in previous experiments 
. The primary beam of 1-A GeV 208 Pb was delivered 
by the heavy-ion synchrotron SIS at GSI, Darmstadt. 
The proton target was composed of 87.3 mg/cm 2 liquid 
hydrogen Jl(| enclosed between thin titanium foils of a to- 
tal thickness of 36 mg/cm 2 . The primary-beam intensity 
was continuously monitored by a beam-intensity monitor 
(SEETRAM) based on secondary-electron emission. In 
order to subtract the contribution of the target windows 
from the measured reaction rate, measurements were re- 
peated with the empty target. Heavy residues produced 
in the target were all strongly forward focused due to the 
inverse reaction kinematics. They were identified using 
the FRagment Separator (FRS) 

The FRS is a two-stage magnetic spectrometer with a 
dispersive intermediate image plane (S2) and an achro- 
matic final image plane (S4) with momentum acceptance 
of 3 % and angular acceptance of 14.4 mrad around the 
beam axis. Two position-sensitive plastic scintillators 
placed at S2 and S4, respectively, provided the magnetic- 
rigidity (Bp) and time-of- flight measurements, which al- 
lowed to determine the mass-over-charge ratio of the par- 
ticles. In the analysis, totally stripped residues were con- 
sidered only. In the case of residues with the highest nu- 
clear charges (above 65Tb) an achromatic degrader (5.3 
g/cm 2 to 5.9 g/cm 2 of aluminum) was placed at S2 to ob- 
tain a better Z resolution. The elements below terbium 
were identified from an energy-loss measurement in an 
ionisation chamber (MUSIC). The velocity of the identi- 
fied residue was determined at S2 from the Bp value and 
transformed into the frame of the beam in the middle 
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of the target taking into account the appropriate energy 
loss. About 100 different values of the magnetic field 
were used in steps of about 2 % in order to cover all 
the produced residues and to construct the full velocity 
distribution of each residue. 
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FIG. 1. Main correction factors applied to the data. Left 
panel: factors due to secondary reactions in all layers of mat- 
ter of the FRS, including the degrader, f 3ec (dashed line), and 
secondary reactions and beam attenuation inside the hydro- 
gen target, ftar (solid line), versus fragment mass number. 
Right panel: factor corresponding to not fully stripped ions, 
f q , versus fragment atomic number. The gray area around 
each line denotes the associated systematic error. 

The measured counting rates N, attributed to a spe- 
cific isotope, were normalized to the number of projec- 
tiles N p recorded with the beam-intensity monitor and to 
the number of target atoms per area n t . Then, the pro- 
duction cross-sections a pro d are calculated by applying 
several corrections according to the following equation: 



N 



? prod 



n t • N p 



fr ' fe ftr ' fq ' ftar ' fs 



(i) 



The used correction factors arise from: the measured 
dead-time of the data-acquisition system (/ T ), the effi- 
ciency of the detection system (/ e ), the loss of fragments 
due to the fragment-separator transmission (ftr), the loss 
of fragments due to incompletely stripped ions (fq), the 
influence of beam attenuation and secondary reactions in 
the liquid- hydrogen target (ftar), and in the other layers 
of matter inside the fragment separator (/ sec ). The f T , 
f £ and ftr correction factors are directly deduced from 
the experiment with high precision. The dead time was 
measured by the acquisition system and kept below 30 
%. The efficiency of all detectors was estimated directly 
from the obtained data to be higher than 98 %. Since 
the full velocity distribution is constructed for each iso- 
tope from the data of different field settings, transmission 
losses are negligeable in the present experiment. 

The corrections due to incompletely stripped ions (f q ) 
and secondary reactions (ftar and f sec ) depend on the 
fragment type. They are displayed with their associated 
uncertainties in Fig. [l]. f q represents the counting loss 
due to incompletely stripped ions. It is significant only 
in case of the fragments with a high nuclear charge. For 
several isotopes the ratio between fully and incompletely 



stripped ions was determined. It made possible to es- 
timate the loss in counting rate due to the fraction of 
incompletely stripped ions for all isotopes. f sec corre- 
sponds to the loss of residues through secondary reac- 
tions in the thick aluminum degrader and other layers 
of matter in the beam line. It was calculated using two 
different formulas for the total reaction cross-sections, 
developed by Karol (l2) and Benesh [i3fl . The results 
agreed within 5 %. f sec varies from 2 to 1.8 with de- 
creasing mass. However for Gd, Eu, Sm and Pm whose 
cross-sections were collected without the degrader, it is 
not higher than 1.13. Secondary reactions inside the hy- 
drogen target also lead to a reduction of the counting 
rates of the heaviest isotopes, but on the other hand pro- 
duce more lighter isotopes. The corresponding correction 
factor, ftar, was estimated from reaction rates obtained 
in the present experiment using a deconvolution method. 
The total reaction cross-section formula of Benesh et al. 
|]l3| utilized in these calculations was adjusted to the ex- 
perimental data from the Barashenkow compilation [pd[ . 

All uncertainties of the used corrections lead to a fi- 
nal systematic error of 9 % to 23 % for Z from 82 to 
61. The measured production cross-sections of the spal- 
lation residues in the reaction of 1A GeV 208 Pb with 
protons are plotted as isotopic distributions in Fig. |^. 
Most of the presented distributions exhibit a Gaussian- 
like shape where the neutron-proton evaporation com- 
petition determines the position of the maximum. The 
most significant deviations from this shape occur for the 
neutron-rich fragments with masses close to that of the 
projectile. In the case of these residues, one and a few 
neutron-removal channels from low excited nuclei created 
mainly in peripheral collisions are responsible for the in- 
creased production cross-sections. Most of the produced 
isotopes populate a corridor, between the valley of sta- 
bility and the proton drip line due to the fact that the 
excited heavy prefragment evaporates mainly neutrons. 

TABLE I. Comparison of the cross-sections from the 
present work (up us) with those obtained by gamma spec- 
troscopy (ctrc) by Gloris et al. The orc values are given 
with total error. The ofrs cross-sections are with statistical 
and total (in parentheses) error. In the fourth column the 
systematic uncertainty for each ofrs is given. 



Isotope orc (mb) 



ofrs ( m b) Systematic error (%) 





22.3±6.1 17.0±0.4(1.6) 


9 


196 Au 


3.88±0.47 4.0±0.1(0.4) 


9 


194 Au 


6.85±0.92 6.3±0.2(0.6) 


9 


148 Eu 


0.104±0.04 0.075±0.005(0.010) 


12 


144 Pm 


0.068±0.013 0.036±0.003(0.006) 
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In Fig. the cross-sections obtained by gamma spec- 
troscopy B arc also shown. To compare with our data, 
we have chosen only isotopes shielded by long-lived or 
stable precursor in the decay chain. A more detailed com- 
parison is presented in Table |. In the case of 196 Au, the 
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FIG. 2. Isotopic production cross-sections of elements between Z—82 and 61, in the reaction of 1-A GeV 208 Pb on hydrogen, 
versus neutron number. Stable (resp. radioactive) isotopes are marked by open (resp. full) triangles. Gamma-spectroscopy 
data regarding shielded isotopes from J6| are plotted as open circles. The solid, dashed and dotted curves were calculated with 
the Cugnon- Schmidt | po|j2l| , Bertini jL6|- Dresner and Isabel jlTj-Dresner models, respectively. 



cross-section is the sum of the production of the ground 
and the isomeric states. The data agree within their error 
bars, except for the isotope with the lowest cross-section, 

144p m 

Spallation reactions are generally modeled as a two- 
step process. In the first step, the nucleon-nucleon col- 
lisions inside the nucleus induce the loss of a few nucle- 
ons and lead to the formation of an excited prefragment. 
This process can be described by the intranuclear cas- 
cade model (INC) sometimes including a pre-equilibrium 
emission. In the second step, the prefragment deexcites 
by evaporation of light particles or by fission. Calcu- 
lations performed with different INC plus evaporation- 
fission models are shown together with our results in 
Figs. H and [| The first two calculations were done with 
the commonly used LAHET Code System (version 2.7 
with default options) from Los Alamos using either 
the Bertini Jl6| plus pre-equilibrium (dashed line) or Is- 
abel |l7j (dotted line) INC models followed by the Dres- 
ner evaporation-fission model . The shapes of the 
isotopic distributions obtained with both INC models are 
very similar and differ significantly from the experimental 
ones: they are shifted with respect to the experimental 
ones towards the neutron-rich side. This can be ascribed 



to the fact that the prediction of the neutron-proton 
evaporation competition in the Dresner code is not sat- 
isfying. The magnitude of the measured and calculated 
cross-sections is also quite different, especially in the case 
of the lighter elements. This effect is better visible on 
the mass distribution (Fig. ^|, upper panel) and more 
marked with the Bertini model which overpredicts largely 
the production of light isotopes. This discrepancy of the 
Bertini model is due to a distribution of excitation en- 
ergies (E*) of the prefragments extending to too high 
values, which results in evaporating more particles and 
finally producing lighter nuclides. This problem of a too 
high E* at the end of the Bertini INC model was already 
noticed in a comparison with neutron double-differential 
cross-section measurements Q although, here, the use 
of the pre-equilibrium option has led to somewhat smaller 
E*. On the other hand, in a region very close to the pro- 
jectile mass, both Bertini and Isabel calculations are in 
good agreement with the data. 

The third calculation (solid line in Fig. || and |J) was 
performed with the version INCL3 of the Cugnon model 
[20 1 combined with a model elaborated by Schmidt et al. 
[ 21 1 . This calculation reproduces much better than the 
former ones the shape of the experimental isotopic dis- 
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tributions. This comes mainly from a better description 
of the neutron-proton competition in the Schmidt than 
in the Dresner evaporation model, since the E* distribu- 
tion at the end of the INC stage is similar in the Isabel 
and Cugnon models (except for very small E*). For ele- 
ments from 760s to 79AU the code predicts shoulders on 
the neutron-rich side of the isotopic distributions. We 
attribute these to the statistical treatment of the Pauli 
blocking in the Cugnon model which improves signifi- 
cantly the excitation-energy distribution p^ j in general 
but also leads to a few prefragments with unrealistically 
low excitation energies. This problem is already partly 
cured in INCL3 compared to the previous version [pO| . 
The magnitude of the cross-sections is not always repro- 
duced, the calculation under-predicting the production 
of the light isotopes. Besides, the main defect of this cal- 
culation is the underproduction of isotopes very close to 
the projectile, which represent an important part of the 
total cross-section. This is ascribed to the sharp surface 
approximation in the Cugnon model which leads to a bad 
description of the most peripheral reactions. These de- 
fects result in a poor prediction of the mass distribution. 



about 8 % to 30 % for A from 140 to 208. Calcula- 
tions performed with the Cugnon-Schmidt (solid line), 
Bertini-Dresner (dashed line) and Isabel-Dresner (dotted 
line) codes are also shown, together with an empirical pa- 
rameterization (dash-dotted line) describing the average 
longitudinal momentum transfer distributions derived by 
Morrisscy from a large compilation of experimental data 
p3j. The Cugnon-Schmidt code predicts recoil energies 
up to 35 % higher than the experimental ones while the 
three other calculations underestimate the data for large 
mass losses. 

In conclusion, the fragment-separator facility at GSI 
has been used to determine, for the first time, the produc- 
tion cross-sections and momentum distributions of 446 
isotopes from spallation reactions of 1-A GeV 208 Pb with 
protons. The results agree with most of the few cross- 
sections previously measured by gamma-spectroscopy. 
Calculations using different models have been performed. 
Although none of them provide a detailed description of 
the data, the new Cugnon-Schmidt code gives clear im- 
provements. 
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FIG. 3. Mass distribution (upper panel) and recoil ki- 
netic energy (bottom panel) of the residues produced in 1-A 
GeV 208 Pb on hydrogen reactions (triangles) versus mass 
number, compared with the Cugnon-Schmidt (solid line), 
Bertini-Dresner (dashed line) and Isabel-Dresner (dotted line) 
models. The dash-dotted line shows the recoil kinetic energies 
expected from the Morrissey systematics p3| . 

The velocity distribution of each residue was also de- 
termined, from which it was possible to infer informa- 
tion about the recoil kinetic energy in the projectile sys- 
tem. In the bottom part of Fig. |^, the average recoil 
kinetic energy of the fragments is shown as a function of 
their mass number. The systematic uncertainty of the 
obtained values (not shown in the picture) varies from 
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